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Abstract 

Recent results from the NuTeV experiment at Fermilab and the deviation of the Z invisible 
width, measured at LEP/SLC, from its Standard Model (SM) prediction suggest the suppression 



Ph, of neutrino-Z couplings. Such suppressions occur naturally in models which mix the neutrinos 

Oh' 

with heavy gauge singlet states. We postulate a universal suppression of the Zvv couplings by a 



factor of (1 — e) and perform a fit to the Z-pole and NuTeV observables with e and the oblique 



^ . correction parameters S and T. Compared to a fit with S and T only, inclusion of e leads to a 

h : 

dramatic improvement in the quality of the fit. The values of S and T preferred by the fit can be 
obtained within the SM by a simple increase in the Higgs boson mass. However, if the W mass 
is also included in the fit, a non-zero U parameter becomes necessary which cannot be supplied 
within the SM. The preferred value of e suggests that the seesaw mechanism may not be the reason 
why neutrinos are so light. 
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I. INTRODUCTION 

The Neutrinos at the Tevatron (NuTeV) experiment 1] at Fermilab has measured the 
ratios of neutral to charged current events in muon (anti)neutrino - nucleon scattering: 

_ a(u,N -> v,X) _ a , 

cr(z/ M X -> /x + X) r 

where 

= cjyV -> ^+X) 1 
' a(^X -> //-X) ~ 2 ' 

and has determined the parameters g\ and g^. |2j to be 

g 2 L = 0.30005 ±0.00137, 

g\ = 0.03076 ±0.00110 . (3) 

The Standard Model (SM) predictions of these parameters based on aglobal fit to non- 
NuTeV data, cited as [g 2 L ] SM = 0.3042 and [g 2 R } SM = 0.0301 in Ref. |l|, differ from the 
NuTeV result by 3a in g\. This disagreement between NuTeV and the SM (as determined 
by non-NuTeV data) is sometimes referred to as the NuTeV 'anomaly' Q]. 



(2) 



Various suggestions have been forwarded as the cause of this discrepancy 



DQDQQ 



Hj. These include theoretical uncertainties due to QCD effects, tree level and loop effects 
resulting from new physics, and nuclear physics effects. We refer the reader to Ref. [3J for 
a comprehensive review. In this paper, we investigate one explanation which is particularly 
attractive in its simplicity and economy. 

Note that the NuTeV value of g\ is smaller than the SM prediction. This is because the 
ratios R u and R p were both smaller than expected, i.e. the neutral current events were not as 
numerous as predicted by the SM when compared to the charged current events. In addition, 
the invisible width of the Z measured at the CERN/SLAC e + e _ colliders LEP/SLC is also 



known to be 2a below the SM prediction jl2j . Both observations suggest that the couplings 
of the neutrinos to the Z boson are suppressed with respect to the SM. 

Suppression of the Zvv couplings occurs naturally in models which mix the neutrinos 
with heavy gauge singlet states 0,0, [3, 14, Q 13- For instance, if the SU(2)l active 



neutrino vl is a linear combination of two mass eigenstates with mixing angle 9, 

U L = (COS flight + (Sin 6>)f heavy , (4) 



then the Zvv coupling will be suppressed by a factor of cos 2 9 if the heavy state is too massive 
to be created in the interaction. Note that the Wtv coupling will also be suppressed by a 
factor of cos 8. 

In general, if the Zvv coupling of a particular neutrino flavor is suppressed by a factor of 
(1 — e), then the W£v coupling of the same flavor will be suppressed by a factor of (l — e/2). 
For simplicity, and to preserve lepton universality, assume that the suppression parameter 
e is common to all three generations. The theoretical values of R u and Rp are reduced by 
a factor of (1 — e), since their numerators are suppressed over their denominators, and the 
invisible width of the Z is reduced by a factor of (1 — 2e). At first glance, then, it seems 
that neutrino mixing may explain both the NuTeV and invisible width discrepancies. 

However, more careful consideration appears to veto this possibility. One of the inputs 
used in calculating the SM predictions is the Fermi constant Gf, which is extracted from 
he muon decay constant G^. Suppression of the W£v couplings leads to the correction 



Gf — G^{1 + e) , (5) 

which will affect all SM predictions. Ref. [j] concludes that a value of e large enough to 
explain the NuTeV anomaly would affect the predictions of the SM (with a conventional 
light Higgs scalar) to an extent that the excellent agreement between the SM and Z-pole 
observables would be lost. 

Notice, though, that the Fermi constant Gf always appears multiplied by the p-parameter 
in neutral current amplitudes. Thus, the predictions for Z-pole observables will be undis- 
turbed if a shift in Gp is compensated by a shift in p. Such shifts can arise via oblique 
corrections due to new physics 

In the following, we perform a fit to the Z-pole and NuTeV data with the oblique correc- 
tion parameters 5* and T as well as the Zvv suppression parameter e. Oblique corrections 
alone cannot explain the NuTeV and invisible width discrepancies. However, fits involving 
the parameter e show excellent agreement with both the Z-pole and NuTeV data. 

The question then, is: "What physics would account for the values of 5* and T that allow 
for a non-zero e?" A simple (though perhaps not unique) solution is the SM itself with a 
heavy Higgs boson. Indeed, we show that if the Higgs boson mass is allowed to be large, 
then the Z-pole and NuTeV data can be fit by e alone. 

The preferred value of e suggests large mixing angles between the active neutrinos and 
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TABLE I: The correlation matrix of the Z-lineshape parameters from LEP. 

the heavy sterile states. This may rule out the seesaw mechanism Q| as the explanation of 
the smallness of neutrino masses. 

If the W mass is included in the fit, we must also introduce the oblique correction param- 
eter U. The preferred value of U is virtually independent of the reference Higgs boson mass 
and cannot be supplied within the SM. We discuss what type of new physics may generate 
the required value of U. 



II. THE DATA 



We begin by listing the data which we will use in our analysis. 
The Z lineshape parameters from LEP/SLC assuming lepton universality are (Ref. 
page 8) 

M z = 91.1875 ± 0.0021 GeV, 

T z = 2.4952 ±0.0023 GeV, 

al = 41.540 ± 0.037 nb , 

R° e = 20.767 ±0.025 , 
A°/ B = 0.0171 ± 0.0010 , (6) 



with the correlation matrix shown in Table |l] Of these parameters, Mz is used as input 
to calculate the SM predictions, and is used to fix the QCD coupling constant ot s {Mz)- 
The remaining three are equivalent to (Ref. [3], pages 8, 9, and 146) 

ri ep t = 83.984 ± 0.086 MeV , 
Tinv/riept = 5.942 ± 0.016 , 
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0.23099 ±0.00053 



(7) 



There is a correlation of 0.17 between r^pt and r inv /Ti ep t w' 
gible. The SM prediction for the Z invisible width is (Ref. 



lile other correlations are negli- 
rj, page 9) 



Tinv/Tlept 



5.9736 ±0.0036 . 



(8) 



As mentioned in the introduction, this is 2a above the experimental value. 

The remaining observables from LEP and SLC, i.e. various asymmetries and ratios of 
partial widths, can be interpreted as measurements of sin 2 9^ in the absence of vertex 
corrections from new physics The average of the values obtained from r polarization, 
^4lr, -^fb) ^fb> an d (Qfb), and that obtained from A F ' B above, is (Ref. 12], page 146) 



sin 3 05* 



d l T = 0.23148 ± 0.00017 . 



(9) 



It should be noted that the agreement among the values of sin 2 B 1 ^ 1 obtained from the six 
observables is not good. The % 2 /d.o.f. associated with the average is 10.2/5. However, since 
neither oblique corrections nor neutrino mixing has any effect on the quality of the agree- 
ment, we will just use the average value, Eq. (jHJ), as representative of these measurements. 
The observables not included in this average, such as and R c , depend only weakly on 
sin 2 and carry little statistical significance [2(| . Indeed, we have also performed an 
analysis with all the asymmetries and heavy flavor data, together with their correlations, 
included separately and have confirmed that doing so does not change the conclusions of 
this paper. 

For g\ and g R , we take the values in the 1998 Particle Data Book [2^, which are based 
on all experiments prior to NuTeV, 



g \ = 0.3009 ±0.0028 
g 2 R = 0.0328 ± 0.0030 



and average them with the NuTeV values, Eq. 0. We obtain, 



(10) 



g \ = 0.3002 ±0.0012 
g\ = 0.0310 ±0.0010 



(11) 



The correlation is negligible. 
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Observable SM prediction Measured Value 



^lept 


83.998 MeV 


83.984 ± 0.086 MeV 
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0.23148 ±0.00017 




0.3037 


0.3002 ± 0.0012 


A 


0.0304 


0.0310 ± 0.0010 


M w 


80.375 


80.449 ± 0.034 GeV 



TABLE II: The observables used in this analysis. The SM predictions are for inputs of M top = 
174.3 GeV, M ffig gs = 115 GeV, a s {M z ) = 0.119, and Aa£jj = 0.02761. They differ from the SM 
predictions cited in the text since the global fit prefers a slightly different set of input values. 



We will also use 
and LEP2 is (Ref. 



;he W mass later in the analysis. The world average from pp experiments 
3, Page 150) 

M w = 80.449 ±0.034 GeV . 



(12) 

The inputs to ZFITTER are the Z mass from Ref. 12 1, listed in Eq. (|Hjl. the top mass from 



The SM predictions are calculated with ZFITTER v6.36 



Ref. 



24 



22 1 and the formulae in Ref. 



M top = 174.3 ±5.1 GeV 



(13) 

and the hadronic contribution to the running of the QED coupling constant from Ref. j^, 

Aa^ d (M z ) = 0.02761 ± 0.00036 . (14) 

The QCD coupling constant is chosen to be 

a s {M z ) = 0.119 . (15) 

All our observables are either purely leptonic, hence QCD corrections only enter at the two 
loop level, or ratios of cross sections in which the leading order QCD corrections cancel. Con- 
sequently, the results of the fits depend only very weakly on this choice. Finally, the Higgs 
mass is varied from 115 GeV, the lower bound from direct searches j^fj, up to 1000 GeV. 
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III. THE CORRECTIONS 



We now consider the effect of the oblique correction parameters S, T, U, and the Zvv 
coupling suppression parameter e on the electroweak observables. The dependence of the 
observables on the oblique correction parameters was obtained in Ref. j^| so we will not 
reproduce it here. To obtain the dependence on e, we note that the tree-level relation 
between sin 2 W and the input parameters a, Gp, and M z is given by 

1/2 



sin 20 



Air a 



y/2G F M 2 z 

Thus, a shift in Gp, Eq. (jSJ), will shift sin 2 9 W (be it sin 2 9^ shc11 ^ 

5Gp 



(16) 



or sin 



l T) by 
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Combining with the oblique corrections, we obtain 



2 2 
S C 



e . 



(17) 



5M W 
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2 121 



4(c 2 - s 
aS 



(aT + e) 



+ 



c — s 
c 2 

lie 1 - s 2 ) 



aT + ^e) + 



aU 



(18) 



A(c 2 -s z ) ' 2(c z -s z ) V ' c z ~) ' 8s 2 
The Z widths are proportional to the product pGpM z (^tJ and it is easy to see that 

pGpM z 1 ' 

From Eqs. (|18|) and ()19jl we see that all Z-pole observables, except for the invisible width, 
depend on e only through the combination (aT + e) . Consequently, any shift in the Z-pole 
observables due to a non-zero e could be compensated by a shift in the T parameter with 
the result that the observables remain unchanged. Finally, the Z invisible width must be 
corrected by a factor of (1 — 2e) and the NuTeV parameters g\ and g\ by a factor of (1 — e). 
Numerically, the observables are corrected as follows: 

Tiopt 



ITlept]sM 
Finv/Tlept 

[Finv/riept]sM 
^lept 



sin 2 ^ 



[sin 2 ^f]sM 

gl 



[fl'IlsM 



1 - 0.0021 S + 0.0093 T + 1.2 e , 
1 + 0.0021 S - 0.0015 T - 2.2 e , 
1 + 0.016 S- 0.011 T- 1.4 £ , 
1 - 0.00905 + 0.022T-0.17e , 
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ff— = 1 + 0.031 S- 0.0067T-3.9e, 



[S'rJsm 

' Un ' - = 1 — 0.0036 S + 0.0056 T 



[Mw\SM 

+0.0042 £7 + 0.22 e . (20) 

Here, [*]sm is the usual SM prediction of the observable * using as input. Note that 
the coefficient of e in g\ is small. This means that the dependence of g\ on e due to the 
suppression of the Zpu and Wtv couplings is almost canceled by the dependence of the 
SM prediction through G^. Ironically, although the original motivation for introducing the 
parameter e was to fit g\, it turns out that g\ is the least sensitive to e of all the observables 
considered. 



IV. THE FITS 

Using the data and formulae from the previous sections, we perform several fits. Initially, 
we exclude the W mass and perform a fit to the remaining five observables only since 
including the W mass serves only to determine the U parameter and does not affect the 
values of the other fit parameters. 

A. Oblique Corrections Only 

First, we perform a fit with only the oblique correction parameters S and T. Using 
M top = 174.3 GeV, MHi ggs = 115 GeV as the reference SM and fitting the expressions in 
Eq. (|2T)J) to the three Z-pole observables and the two NuTeV observables we obtain 

S = -0.09 ±0.10 , 

T = -0.13 ±0.12 , (21) 

with a correlation of 0.89. The quality of the fit is unimpressive: x 2 = 11.3 for 5 — 2 = 3 
degrees of freedom. The preferred region on the S-T plane is shown in Fig. 1. As is evident 
from the figure, there is no region where the la bands for Fi ept , sin 2 6^*, and g\ overlap. 

Fig. 1 also shows the preferred region if S and T are fit to the Z-pole observables only. 
Including the NuTeV data shifts both central values of S and T to the negative side. Though 
the g\ band tries to pull S and T into the fourth quadrant (which corresponds to the SM 
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1.00 
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0.56 


1.00 
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-0.20 


-0.73 1.00 
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0.18 


-0.64 0.58 1.00 



TABLE III: The correlations among the oblique correction parameters and e. The correlations 
among S, T, and e are unaffected by whether M\y and the U parameter are included in the fit. 

with a larger Higgs mass), the narrow sin 2 9^ band forces S and T to move along it into 
the third quadrant. 

In conclusion, oblique corrections alone cannot explain the NuTeV anomaly and/or the 
invisible width discrepancy. 



B. With Neutrino Mixing 

Next, we perform a fit with S, T, and e. The reference SM is M top = 174.3 GeV, 
^Higgs = 115 GeV as before. The result is 

S = -0.03 ±0.10 , 
T = -0.44 ±0.15 , 

e = 0.0030 ±0.0010 . (22) 

The quality of the fit is improved dramatically to x 2 = 1.17 for 5 — 3 = 2 degrees of 
freedom. The correlations among the parameters are shown in Table IIHI The preferred 
regions in the S-T, S-e, and T-e planes are shown in Figs. 2 through 4. The central 
values have shifted somewhat from our preliminary report in Ref. 28| . This is partly due 
to our use of updated NuTeV numbers from the third reference of and partly due to 
our use of Aa^ A (M z ) = 0.02761 from Ref. H instead of the ZFITTER default value of 
A«[ 5 i(M z ) = 0.02804. 

Thus, by including both oblique corrections and e, we obtain an excellent fit to both the 
Z-pole and NuTeV data. But what kind of 'new physics' would provide such values of S 
and T? While the limits on S permit it to have either sign, T is constrained to be negative 
by 3a. Few models of new physics are available which predict a negative T 



Recall that the effect of a SM Higgs heavier than our reference value (here chosen to 
be 115 GeV) is manifested as shifts in the oblique correction parameters. The approximate 
expressions for these shifts are [3^ 

Sffiggs ~ ^ln(M Higgs /Mgf ggs ) , 
Tmggs « -^ln(M Hig g S /M- f ggs ) , 

f/uiggs « . (23) 

Notice that Sm ggs is positive while T H ig gs is negative for M Higgs > Mfj e iggs . Thus increasing the 
Higgs mass will have the desired effect of providing a negative T. Can T be made negative 
enough without making S too large? The answer is provided in Fig. 5, in which we show 
a blow up of the central region of Fig. 2. The SM points fall comfortably within the 90% 
confidence contour when the Higgs mass is even moderately large. 



C. Neutrino Mixing Only 

To check that the SM itself is indeed compatible with the data (excluding the W mass) 
we perform a fit with only e as the fit parameter, and plot the dependence of x 2 on the 
Higgs mass used to define the reference SM. The result is shown is Fig. 6. We also show the 
^Higgs dependence of the la limits on e in Fig. 7. The graphs demonstrate that the data 
are well fit by a SM with a heavy Higgs and a nonzero value of e. 



D. The W Mass 



We have deliberately excluded the W mass from our analysis since it cannot be fit with 
e alone, or by S, T, and e. Eq. (fTHj) suggests that the single parameter fit with e, which 
yields a positive value of e, may increase the W mass prediction towards the experimental 
value 4]. However, the coefficient of e, Eq. (|20J) . is too small for it to significantly mitigate 
the discrepancy. A three parameter fit with S, T, and e also fails to close the gap. In that 
case, note that the combination of T and e upon which My/ depends is 

aT + ^ = (oT + e) - e ■ (24) 

Since (aT + e) is pinned by the Z-pole observables, a non-zero e would actually lead to a 
decrease in the W mass prediction. Therefore, the U parameter is necessary to fit My/. The 
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result of the four parameter fit to the reference SM of M top = 174.3 GeV, Mmggs = 115 GeV 
is 

U = 0.62 ± 0.16 , (25) 

and the correlations with the other parameters are shown in Table 11111 The dependence of 
these limits on the reference Higgs mass Mmggs is weak, as implied by Eq. (|2*3|) and shown 
in Fig. 8. 



V. CONCLUSIONS AND DISCUSSION 



We conclude that neutrino mixing together with oblique corrections can reconcile the 
Z-pole and NuTeV data. The simplest model which provides the necessary values of S and 
T is the SM itself with a large Higgs boson mass. 

We emphasize that there is nothing exotic about this solution. Neutrinos are known 
to mix from Super-Kamiokande 

0, SNO K2K H, 

and other neutrino experiments. 
Neutrino mixing may have the added bonus of contributing negatively to both S and T . 

The possibility that the Higgs boson is heavy has also been under scrutiny [34I I35I since 
it has not been found in the ~ 80 GeV range preferred by the SM global fit (Ref . [12j , page 
154). If the Higgs boson is indeed heavy as suggested by Fig. 6, it may favor a possible 
dynamical mechanism for electroweak symmetry breaking while providing a challenge for 
super symmetry j^. 

The limits on the suppression parameter, 

e = 0.0030 ±0.0010 , (26) 

imply 

6 = 0.055 ±0.010 , (27) 

as the mixing angle if the suppression is due to mixing with a single heavy sterile state. The 
heavy state into which the neutrino mixes must be at least as massive as the Z so that the 
Z does not decay into it. A naive seesaw model 18] does not permit such a large mixing 
angle since the angle is related to the ratio of masses. However, it can be shown that the 
required pattern of mixings and masses can be arranged when there exist inter-generational 
mixings [15]. We will present explicit examples in Ref. 
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The W mass needs the U parameter, so we cannot do without other new physics entirely. 
The Higgs mass dependence of the oblique correction parameters are shown in Fig. 8. What 
kind of new physics will be compatible with this? It has to predict a small T while predicting 
a large U. One possibility is that the U parameter can be enhanced by the formation of 
bound states at new particle thresholds. If one expresses T and U as dispersion integrals 
over spectral functions, one finds 

r°° ds 



— [ImII±(s) - Imllo(s)] , 

-thrcs ^ 

U oc / — [lmn±0) - ImlloOO] , (28) 



where we have used the notation of Ref. 



381 ] . Due to the extra negative power of s in the 



integrand of U, it is more sensitive to the enhancement of the threshold than T. Indeed, it 
has been shown in Ref. j3] that threshold effects do not enhance the T parameter. This 
could, again, mean that a theory that exhibits dynamical electroweak symmetry breaking is 
favored. This possibility will be investigated in a subsequent paper. 
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FIG. 1: The fit to the data with only S and T. The bands associated with each observable show the 
lcr limits. The shaded ellipses show the 68% and 90% confidence contours. The unshaded ellipses 
partially hidden behind the shaded ones show the contours when only the LEP/SLC data is used. 
The origin is the reference SM with M top = 174.3 GeV and -Mmggs = 115 GeV. The curved arrow 
attached to the origin indicates the path along which the SM point will move when the Higgs mass 
is increased from 115 GeV to 1 TeV. The curved lines above and below it indicate the uncertainty 
in the SM point's position from that of the top mass: M top = 174.3 ±5.1 GeV [24]. 




FIG. 2: The 68% and 90% confidence contours on the S-T plane when the data is fit with S, T, 
and e. The SM points are shown as in Fig. 1. 
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FIG. 3: The 68% and 90% confidence contours on the S-e plane. The arrow attached to the origin 
indicates the path along which the SM point will move when Mmggs is increased from 115 GeV to 
1 TeV. The dependence of the SM point on Mt op is not shown. 




FIG. 4: The 68% and 90% confidence contours on the T-s plane. The arrow attached to the origin 
indicates the path along which the SM point will move when MHi ggs is increased from 115 GeV to 
1 TeV. The dependence of the SM point on M top is not shown. 
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FIG. 5: Blow up of the central region of Fig. 2. The Mm gg and M top dependence of the SM point 
is shown as in Fig 2, with the Higgs mass indicated in units of GeV. The horizontal arrow going 
through the origin indicates the uncertainty in the position of the SM point due to the uncertainty 
in Aa^ d (Mz), with the arrow pointing in the direction of increasing Aa^ d (Mz) 
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FIG. 6: The dependence of x 2 /d.o.f. on MHi ggs when only e is used to fit the Z-pole and NuTeV 
data. The W mass is not included in the fit. The Higgs mass in the shaded region is excluded by 
direct searches. 
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